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Abstract. We present a simple and accessible method wiiicii uses contour 
integration methods to derive formulae for functional determinants. To make 
the presentation as clear as possible we illustrate the general ideas using the 
Laplacian with Dirichlet boundary conditions on the interval. Afterwards, we 
indicate how more general operators as well as general boundary conditions 
can be covered. 



1. Introduction 

In the Euclidean path-integral formalism, physical properties of systems are con- 
veniently described by means of path-integral functionals. The leading order con- 
tribution to path-integrals often involves the integration of the exponential of a 
quadratic form. This integration can be carried out exactly and yields a functional 
determinant U 12] ■ 

In recent years, a contour integration method for the evaluation of determinants 
has been developed PJI^IFj. Although the applicability of the method in higher di- 
mension seems to be restricted to highly symmetric configurations, it turns out that 
in one dimension quite generic situations can be considered 6 . In particular, the 
case with zero modes present can be dealt with very elegantly and the determinants 
with the zero mode extracted can be found in closed form. 

We will first introduce the contour integral methods using the simple example 
of the Laplace operator with Dirichlet boundary conditions on the interval. There 
will be no zero modes present and the example serves to introduce and explain 
some details of the method. Based on the example, more general second order 
operators with Dirichlet boundary conditions will be considered. The results found 
can be compared with known answers|3|S]. Afterwards, we explain the additional 
complications when zero modes are present. The conclusions summarize the main 
results and will indicate further possible generalizations of the approach to Sturm- 
Liouville type operators with arbitrary boundary conditions. 



2. Example for the Laplacian on the interval 

Let us consider first the Laplacian L = on the interval / = [0,1] with 

Dirichlet boundary conditions. The eigenvalues A^j and eigenfunctions Un are de- 
fined by the equations 

(1) - -j^Unix) = Xlu„{x), 

(2) Un{0) = U„(1)=0. 
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Equation ^ is satisfied by 

(3) Un{x) = Ci cos(A„a;) + C2 sin(A„a;). 

The boundary condition u„(0) = forces Ci = 0, and u„(l) =0 forces 

sinA„ =0 <^ A„ = 7rn, ri £ "Z. 

The set of nontrivial Hnearly independent eigenfunctions is given by n e IN. For 
this case where eigenvalues are known explicitly, it is straightforward to write down 
the zeta function associated with L, 

oc 

ri=l ^ ' 

and to calculate the functional determinant in the zeta function schenie|^ llOllll| . 

(4) - In dot L = Cl (0) = -2Cfi(0) hi tt + 2Cr{0) = In tt - ln(27r) = - In 2. 

However, for cases where contains a potential, eigenvalues are mostly not known 
and a simple procedure such as the one above is not available. To make this more 
general case accessible, we will therefore introduce contour integration methods to 
study the above example. The methods are then directly generalizable to cases 
with a potential. 

Starting point for the contour integration method is the implicit eigenvalue equa- 
tion 

(5) sinA„ — 0. 

This equation allows the zeta function to be written as the contour integral 



dk k Insinfc, 
ok 



where the contour 7 is counterclockwise and encloses all eigenvalues on the positive 
real axis. As given, the representation is valid for 5Rs > 1/2. Next, one would 
like to shift the contour towards the imaginary axis. In doing so, we pick up a 
contribution from the origin because sinO = 0. Therefore, a slightly more elegant 
formulation is obtained when starting with 

7 

where the contribution from the origin is avoided because now limfe^o = 1 . So 
now, nothing prevents us from shifting the contour to the imaginary axis and by 
taking due care of the cut of the complex root, which we place along the negative 
real axis, we find 

00 

^ , , siuTTs f „ , „2<! 9 , sinhfc 
(7) ais) = ^jdkk -^In^. 



As the fc — > behavior, and the A: — > 00 behavior, obtained from 

sinhfc 



= ^(1 



fc 2k 
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show, this integral representation is weU defined for 1/2 < 5fts < 1. The restriction 
SRs < 1 comes from the k ^ behavior and 5Rs > 1/2 comes from the k ^ oo 
behavior. Writing for some e > 



Sm TTS 



,, , 9„ 9 , sinhfc sin TTS 

dkk~^' — In — — + 

ok k TT 



dkk-''—\n'^^ 
dk^"" k ' 



we need to construct the analytical continuation of the second integral to s = 0. 
This is achieved by adding and subtracting the leading fc — > oo asymptotics, thus 
writing 



Cl(s) 



sm TTS 



dfcfc-^"— In 
ok 



d , sinhfc 



sm TTS 



sm TTS 



dfcfc"'''*— In 



d , / sinh k 



dkk-^^'—lTl 

ok 



dk \ k 
sinh k 



2ke 



sm TTS 



dkk-^'—ln 
Ok 



sm TTS 



In 



d , / sinh k 



dk \ k 



2ke' 



sm TTS 



2s - 1 



2s 



valid around s = 0. 

The derivative can now be easily carried out and its value at s = is given in 
its simplest form by letting e — > 0. We obtain 



(8) Ci(0) = hm In 



— lim In 



/ sinh k 
/ sinh k 



In 



In 2 



^'J^2ke-' 
In 2, 



In, 



as given in Q). 



3. More general second order operators 



Generalizing from the previous example, it is the aim of this section to develop 
the contour integration approach for operators of the type 

d2 



(9) 



L 



+ V{x), 



again on the interval / = [0,1]. The potential V{x) is assumed to be continuous 
and to start with all eigenvalues will be assumed to be positive, so no zero modes 
are present. With an eye on equation © we let Uk(x) be the unique solution of 

(10) (L-fc2)ufe(x) =0 
satisfying 

(11) ufc(o) = o, 4(0) = 1. 

The condition Ufc(O) = is chosen such that the eigenvalues are fixed by imposing 

(12) Mfe(l)=0. 
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Then, the contour integral representation of the associated zeta function is 

(13) a(.)=2^/*fc-^|ln«.(l). 

7 

The normahzation u'^(O) = 1 is chosen such that when shifting the contour towards 
the imaginary axis no contributions from the origin arise, see ©. As the leading 
k ^ oo asymptotics does not depend on the potential V{x), the subsequent steps 
are precisely as for the example in Section 2. Therefore, a glance at equation (jH)) 
shows that 

(14) Cl(0) = -limln(w,fe(l))-ln2 = -lnMo(l)-ln2 = -ln(2uo(l)), 

where, by definition, uq{x) is the solution of the homogeneous differential equation 
Luo{x) = with uo(0) = 0, Uo(0) = 1. 

4. Presence of zero modes 

For the case when zero modes are present, the formula H14(l can not be applied 
simply because, as a zero mode, uq(x) satisfies uq(1) = 0. Therefore, the contour 
cannot be shifted to the imaginary axis; instead, particular attention is needed when 
approaching the origin in the /c-plane. We first need to determine the behavior of 
Mfe(l) for small k in order to eliminate the pole in the integrand. To do this, we 
note that integrating the left-hand side of 

/ dx uo{x)Luk{x) = k^ dx uo{x)uk{x) = k'^{uo\uk) 
Jq Jo 

by parts gives 

[u'f^{x)ukix) - uJ,(x)uo(a;)]o + / dx Uk{x){JjUQ{x)) = fc^(uo|ufc). 

Jo 

Using the boundary conditions this yields 

(15) Ufc(l) = ^-p— = -k fk. 

Since fk is finite and non-zero as fc ^ 0, we have the desired behavior of u/c(l)- 
The minus sign has been included in the definition of fk for later convenience. It 
is important to note that l|15|l is true for any k. 

We can now modify the discussion of Section 3 to cover the case when a zero 
mode is present, by using the following observations. The function fk, defined by 
(|15|l . vanishes at all values of k for which k'^ is a positive eigenvalue, as Ufe(l) does. 
Only in the case of the zero eigenvalue do we have the situation where /o 7^ 0, but 
uo(l) = 0. Thereforee if we use fk instead of Ufc(l) in ifT^ . and choose the contour 
so as to surround the positive eigenvalues, the contour can be shifted towards the 
origin. However, the notion that Uk{l) may be straightforwardly substituted by fk 
has to be amended slightly, since these functions have different behaviors as fc ^ cx). 
This difference in behavior is accounted for if we replace Ufe(l) not by fk, but by 
(1 — k'^)fk, since then Uik{l) ~ k^fik as A: ^ cxo, while (1 — fc^)/fc remains non-zero 
at fc = 0. 
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These remarks lead us to consider the contour integral 

(16) -^ [ dkk-^^^\n{l-k')h^CLds) + 7^ I dkk-'-^^Hl-k') 

where the contour surrounds all the eigenvalues on the positive k axis. Depending 
on whether the contour encloses the point A: = 1, this might differ from by a 

constant. But given our focus is on Cl{0), this is irrelevant. 

In summary, we are now in a position to proceed precisely as before. In (jl4|l . 
instead of having uo(l) make its appearance, its now /o that comes into play. The 
contributions resulting from the asymptotics at infinity are not influenced at all 
by the existence of a zero mode and therefore remain the same as above in (|14|l . 
Therefore, again referring to this equation, the answer now reads 

(17) Ci(0) - -limln(/fe)-ln2 = -ln/o-ln2 = -ln(2/o). 

For a particular problem specified, namely for any given potential V{x), the de- 
rivative of the zero mode at the endpoint of the interval can be easily determined 
numerically and so can the determinant. 

5. Conclusions 

The calculation provided can be generalized in various respects. First, instead 
of considering iPJ, one can analyze determinants for general Sturm-Liouville type 
operators 

Reformulating the second order differential operator as a first order system we then 
introduce Vk{x) — P{x)u'f.(x). General boundary conditions are imposed as 

where M and N axe 2x2 matrices whose entries characterize the nature of the 
boundary conditions. For this quite generic situation closed formulas for the de- 
terminant can be derived and applied for different examples 12.^. We can also relax 
the condition that the eigenvalues are nonnegative. In order that the associated 
zeta function can be defined in the presence of negative eigenvalues, we place the 
cut of the complex root at some angle different from and from tt. Choosing a 
suitable contour 7 for the representation of the zeta function, only minor changes 
to the procedure explained are necessarv|12|. 
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